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Dedicated to Percy Deift on the occasion of his sixtieth birthday 

Abstract. Unitary random matrix ensembles 

Z~jv(det M) a exp(-A Tr V{M)) dM 

denned on positive definite matrices M, where a > — 1 and V is 
real analytic, have a hard edge at 0. The equilibrium measure as- 
sociated with V typically vanishes like a square root at soft edges 
of the spectrum. For the case that the equilibrium measure van- 
ishes like a square root at 0, we determine the scaling limits of the 
eigenvalue correlation kernel near in the limit when n, N — > oo 
such that n/N — 1 = 0(n~ 2 / 3 ). For each value of a > —1 we find a 
one-parameter family of limiting kernels that we describe in terms 
of the Hastings-McLeod solution of the Painleve II equation with 
parameter a + 1/2. 



I. Introduction and statement of results 

We deal in this paper with the following question. 

Question: Given a unitary random matrix ensemble 

-^(detM) a exp(-NTrV(M))dM, a > -1, (1.1) 

Z n ,N 

defined on positive definite Hermitian matrices M of size nxn, 
where the real analytic potential V is such that the equilibrium 
measure associated with V has a density on [0, oo) that vanishes 
like a square root at 0. What are the scaling limits around of 
the eigenvalue correlation kernel as n, N — + oo, n/N — » 1 ? 

2000 Mathematics Subject Classification. Primary 15A52, 34M55; Secondary 
35Q15. 

The authors are supported by FWO-Flanders project G. 0455. 04, by K.U. Leuven 
research grant OT/04/21, by the Belgian Interuniversity Attraction Pole P06/02, 
and by the European Science Foundation Program MISGAM. The second author 
is also supported by INTAS Research Network 03-51-6637, and by a grant from the 
Ministry of Education and Science of Spain, project code MTM2005-08648-C02-01. 

1 



2 



TOM CLAEYS AND ARNO B.J. KUIJLAARS 



One expects a universality type result for this situation where the 
limiting kernels do not depend on the exact form of V. This is in line 
with known universality results in the usual unitary random matrix 
ensembles of the form 

-^|detM| 2Q exp(-NTrjy(M))dM, a > -1/2 (1.2) 

Z n ,N 

defined on all Hermitian matrices M of size nxn, [2J [31 El [TJ [121 HH] • In 
(11.11) and (II. 2p and throughout the paper we assume that the confining 
potentials V and W are real analytic and satisfy 

V(x) . W{x) 

hm - — — — = +oo, and hm - — — — = +oo. 

x^+oo log(ar + 1J x^±oo log(ar + 1) 

(1.3) 

The random matrix ensemble (11.11) is restricted to positive definite 
matrices, thereby creating a hard edge at 0. This has an influence 
on the local eigenvalue behavior near since eigenvalues are always 
positive. 

It is well-known, see e.g. [HI EE], that the eigenvalue correlation kernel 
for the ensemble (11.11) has the form 

K n>N {x,y) = xhn-^e-^ £pf> (x)pf (y), (1.4) 

where pj^ is the j-th degree orthonormal polynomial with respect to 
the weight x a e~ NV ^ on [0, oo). The limiting mean eigenvalue density 

Vv(z) = lim -K n , )N (x,x) 

n,N — voo yi 
n/N^l 

is the minimizer of the weighted energy 

IyOift) = [[ log- r ip(x)ip(y) dxdy + [ V(x)ip(x) dx (1.5) 

J J \x-y\ J 

taken over probability density functions ip supported on [0,oo), see 
[HI EH]. It does not depend on a. As V is real analytic, the support of 
ipv consists of a finite union of intervals [11] . 

We can then distinguish two situations. The first is that does not 
belong to the support Sy of ipy. Then the hard edge at does not 
lead to a new phenomenon. Generically one finds that the eigenvalue 
correlation kernel (II .4p has the sine kernel 

smnjx-y) 

; r- (1-6) 

tt(x - y) 



WHEN THE SOFT EDGE MEETS THE HARD EDGE 



3 




Figure 1. The density ip Vc for V c (x) = with c 

equal to 0.7 (left figure), 1 (middle figure) and 1.2 (right 
figure). 



as scaling limit in the bulk and the Airy kernel 

Ai(x) Ai'(y) - Ai'(z) M(y) 



;i.7) 



x-y 

as scaling limit at the (soft) edge points of Sy- These are also the 
generic scaling limits for the correlation kernel associated with the ma- 
trix ensemble flO}, see [2j El QUI E2l [21] ■ 

The second situation is that belongs to Sy- In that case, one 
typically finds that ijjy is unbounded at with a square root type 
singularity at 0. Then instead of the Airy kernel we are led to a Bessel 
kernel 

W^i) (L8) 

that describes the local eigenvalue correlations near 0, see [TH1 E21 123] . 

Changing parameters in V we may create a transition from the Airy 
kernel to the Bessel kernel. This happens for example if 

V(x) = V c (x) = ^(x-2) 2 . 

When considered on R this is a shifted and scaled GUE ensemble. For 
c < 1 we have 

Mx) = ^-^/Ac-ix-2) 2 , for x G [2- 2 Vc, 2 + 2^]- 

Z7TC 

These are all semi-circle laws which tend as c — > 1~ to the semi-circle 
law on [0,4]. Then the left soft edge meets with the hard edge at 0. 
For c > 1 we have 

ipv c ( x ) = n ^ 7= ( x + a)Vb - x, for x G [0, b], 
2,-kc^x 



with 



4 4 



-- + -Vl + 3c and 6 = - + -VI + 3c, 
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which has a square-root singularity at 0. See Figure [T] for the density 
tpVc f° r the values c = 0.7, c = 1 and c = 1.2. 

Thus our Question asks about the transitional case. What happens 
when the soft edge meets the hard edge? Note that this transition 
from soft to hard edge is different from the one considered in [4J where 
the authors start from a hard edge situation and let the parameter a 
increase to infinity. Other types of singular edge behaviors in unitary 
random matrix ensembles have been studied in [IB], where a soft edge 
coincides with a spectral singularity, and in [8], where the limiting mean 
eigenvalue density vanishes faster than in the regular case. 

Remark 1.1. Our Question is related to the calculation of the so-called 
Janossy densities for unitary random matrix ensembles 

- — exp(-iV Tr V(M)) dM 

Z n ,N 

in an asymptotic limit. Indeed, it was shown in [5] that the Janossy 
densities for an interval I C R are expressed in terms of the orthogonal 
polynomials for the weight e~ NV ( x ' restricted to the complement of 
I. If tpy vanishes like a square root at a left edge point a, and if 
/ = (— oo,a], then after linear scaling we can reduce this to the case 
a = 0, and then the Janossy densities are expressed in terms of the 
orthogonal polynomial kernel (jl-4p with a = 0. 

Our first result is that there is indeed a one-parameter family of 
limiting kernels that depends on a. 

Theorem 1.2. For every a > — 1, there is a one-parameter family of 
kernels IKa° ft//hard (:r, y; s) depending on s G R, such that the following 
holds. Let V be real analytic on [0, oo) such that ipy vanishes like a 
square root at and such that there are no other singular points in the 
external field V . Then there are positive constants c = c±y and c<iy 
such that 

hm -\j- z K n , N -fr-^J) = K^ d (x, y; s) (1.9) 

n,jv->oo (cnyi 6 \{cn) zli [crip 6 J 

whenever n,N^oo such that 

lim n 2/3 (— -l) = L eR (1.10) 

and s = C2yL. The limit ( 11. 91) holds uniformly for x and y in compact 
subsets of (0, oo). 

See equation (13.31) below for explicit formulas for the constants c\y 
and c 2 y. We also refer to [12] and Section 3 below for the notion and 
the classification of singular points in an external field. 
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In contrast to other works on universality, see e.g. [7J [12], we do 
not have to go through an explicit Riemann-Hilbert steepest descent 
analysis in order to prove Theorem 11.21 Instead we prove Theorem 11.21 
by relating the random matrix ensemble (11.11) to the random matrix 
ensemble (11 ,2p defined on all n x n Hermitian matrices where 

W(x) = ^V(x 2 ) for x e R 

and with a replaced by a ± 1/2. The equilibrium measure ipw for (11.21) 
vanishes quadratically at 0. The universal limiting kernels for this 
situation were found in [7]. They are constructed out of ^-functions 
associated with the Hastings-McLeod solution of the Painleve II equa- 
tion. Using this connection we find the following representation for the 
limiting kernels Ka° ft//hard (x, y; s). 

Theorem 1.3. Let a > —1. Then the kernels Ka° ft//hard (x, y; s) have 
the form 

jj^soft/hard/^, g \ _ fa{ X \ s )9a(y'i s ) ~ fa(y'i s )9a( x i s ) -q\ 

n(x-y) 

where f a and g a are solutions of the system of differential equations 

2q + a/(2x) 2 + (q 2 + r + §)/x\ / f a (x; sf 
-2x-g 2 + r-| -2q-a/(2x) J \g a {x; s) 

(1-12) 

where q = q(s) is the Hastings-McLeod solution of the Painleve II equa- 
tion with parameter a + \ 

q" = sq + 2q 3 -<*-- (1.13) 

and r = r(s) = q'{s). The solutions of (II. 12ft are characterized by the 
asymptotic behavior 




f a (x; s) = x~ l/A cos ( \x^ 2 + sx 1 / 2 - ^(a + 0(x~ z ; - 



;i.i4) 



//„(./■: .s) = -x 1 ^ sin ( ^x 3/2 + sx 112 - ^n(a + \))+ 0{x^ : i 



;i.i5) 



as x — > +oo. 



Remark 1.4. The Hastings-McLeod solution of (11.131) with a = — 1/2 
was introduced in [15]. It is chararacterized by the asymptotic behavior 
q(s) ~ Ai(s) as s ^ +oo where Ai is the usual Airy function. For 
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general a > — 1, a ^ —1/2, the Hastings-McLeod solution of ( 11. 13ft is 
characterized by the asymptotic behavior q(s) ~ (a+ |)/s as s — > +oo 

and g(s) ~ ^J—s/2 as s — ► — oo, see [T31 Chapter 11.7]. It was shown 
in [7] that the Hastings-McLeod solution has no poles on the real line 
and so the system (11.121) is well-defined for every s6R. 
Note that p\ = q 2 + r + | satisfies 

(p') 2 - a 2 
2p 

and P2 = q 2 — r + | satisfies 

„ o 2 , b') 2 -(« + i) 2 



p" = 2p 2 — sp + 



2p 

which are versions of the Painleve XXXIV equation, see [13] . 

The functions f a and g a also satisfy the following linear differential 
equation with respect to s 

d ffa(x;s)\ ( q 1 \ ( f a (x\ 
>9afa;8)) \-x -q) \g a (x;s) 



ds 



(1.16) 



which together with fl 1 . 1 2 [) constitutes the Lax pair for Painleve XXXIV, 
see [131 P a g e 175]. We will not use dl . 16[) in what follows. In [16], a 
limiting kernel is obtained which is built out of solutions to the Lax 
pair for Painleve XXXIV as well. However, the relevant solution to 
Painleve XXXIV in this paper is different from ours, and not related 
to the Hastings-McLeod solution to Painleve II. 

Remark 1.5. The system (j!.12p can also be considered for complex 
values of x. The functions f a and g a have asymptotics in the complex 
x-plane given by 

f a {x; s) = -^ e ^ a+ ^e- i{ l x3/2+sxl/2 \l + 0(aT 1/2 )) (1.17) 

g a ( X - S ) = ^ e §-( a+ §) e -i(§*3/W/ 2 ) (l + (1.18) 

uniformly as x — * oo in the sector e < arga; < 2tt — e for any e > 0. 
The positive real line argx = is an anti-Stokes line for the system 
(I1.12p . where the asymptotic behaviors (I1.17P and (11.181) are not valid. 
As usual one has a two-term asymptotic approximation along the anti- 
Stokes line which leads to (11.141) and (I1.15p . 

Remark 1.6. Besides the eigenvalue correlation kernel, one might be 
interested in the limiting distribution of the (re-scaled) smallest eigen- 
value of the ensemble (II. ip . These distributions are given by Fredholm 
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determinants of the limiting kernels ( II. lip . For general a we do not 
know if these kernels allow Painleve expressions but for a = we can 
relate them to the well-known Tracy- Widom distribution [21] 



F T w(x)=exp(-^ (y-x) 2 q(y)dy^ (1.19) 

where q(y) is the Hastings-McLeod solution of Painleve II. 

Indeed, since V is real analytic on [0, oo) we can also consider V on 
(— 5, oo) for some 5 > 0. For small enough 5 the equilibrium measure 
associated with V on (—5, oo) will not depend on 5 and it will vanish 
like a square root at 0. Then for the ensemble f 1 1.11) with a = and 
defined on Hermitian matrices with eigenvalues > —5 we have that 
the distribution of the smallest eigenvalue tends to the Tracy- Widom 
distribution (11.191) as n — > oo. Here we assume for convenience that 
N = n; otherwise under the assumption (jl.lOp we find a shifted version 
of (11.191) . Moreover, denoting the smallest eigenvalue by Ai we have 
for x > 0, 

P(Ai >x)= P(Ai > x | Ai > 0) = P(Al > X \ (1.20) 

— >0) 



Here P denotes the probability in the ensemble ( 11.11) on positive-definite 
matrices, while P is the probability in the ensemble (II. ip on Hermitian 
matrices with eigenvalues > —5. Hence for an appropriate constant 
c > we have if a = and N = n, 

lim F(cn 2 / 3 X 1 > x) = ^ w( ~^ . (1.21) 

For a^0, the right-hand side of (11.211) will be different, but the scaling 
with n 2//3 in the left-hand side will be the same. 



2. A QUADRATIC TRANSFORMATION 

In this section we relate the correlation kernels (11.41) for the ran- 
dom matrix ensemble (II .ip to the eigenvalue correlation kernel for the 
random matrix ensemble (11.21) with 

W{x) = ^V(x 2 ), xeR (2.1) 
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and parameter a ± 1/2 instead of a. We use K ( ^'^ izl ' W \x, y) to denote 
the correlation kernel. Thus 

j=0 

(2.2) 

where P- is the j-th degree orthonormal polynomial with respect to 
the weight |x| 2a±1 e- w(x) Qn R> 

To distinguish this kernel from the kernel (11.41) that is related to V 
and parameter a we write in this section 

K HtN (x,y) = K%jp{x,y). 

The following result is not new, see e.g. pQ Appendix B] and [TH Exer- 
cises 4.4]. For the convenience of the reader we have included a proof. 

Proposition 2.1. Let V be defined on [0, oo) and let W be defined by 
(12.11) . Assume that (ll.3p is satisfied. Then for a > — 1, we have 

K i$P( x >y) = \( x vY 1a {^S^N^iVx, y/y) + Kf^' W \yfx, -y/y) 

(2.3) 

while for a > 0, we have in addition 

(2.4) 

Proof. In the proof we use the abbreviations 

We also write p^ for the orthonormal polynomial of degree n with 
respect to the weight w(x) = x a e~ NV ^ on [0, oo), and pf^ and 
for the orthonormal polynomials with respect to the weights W\(x) = 
\ x \2a+i e -NW(x) and W2 ( x j = \ x \2a-x e -NW{x) ^ respectively. Note that 

p£^ and Pn are well-defined for a > —1, while is only defined 
for a > 0. 

Since u>i is an even weight on JR., the even degree orthogonal poly- 
nomials are even and the odd degree orthogonal polynomials are odd. 
Therefore there exists a polynomial p£^ of degree n such that 

PST ) (x)=^) (2.5) 



WHEN THE SOFT EDGE MEETS THE HARD EDGE 
For k < n, we then have 

POO 

= / P { l N \x)x 2k \x\ 2a+l e- 2NW ^dx 



oc 
oc 



pW(x 2 )x 2k x 2a+1 e- NV(x2 Ux 



pW(s)s k s a e- NV(s) ds, (2.6) 
where we put s = x 2 . A similar calculation shows that 

/oo n 
(P 2 ( f } (x)) \x\ 2a+1 e- 2NW ^dx 
-oo 



oo 



X 



(p^(s)Ys a e- Nv ^ds, (2.7) 

JO 

Therefore p^ has the orthogonality conditions that characterize the 
orthonormal polynomial p^ and it follows that 

J^(*)=rfV). (2-8) 
If a > 0, we find in a similar way that 

Q { Z\(x)=xpl N \x 2 ). (2.9) 
From (12.21) . we now get for x, y > 0, 

'2n-l 2n-l \ 

E fffc^W" ^ + E 'f^if W5> . 

. 3=0 j=0 / 

(2.10) 

Since p( 2N ^ is even for even j and odd for odd j, we see that the odd 
j-terms in the two sums of ( 12.101) cancel out, while the even j'-terms 
are equal. We then obtain by (12. 8p and (11.41) that 

k L,2n( Vx, Vv) + K^ 2N ( y/E, -y/y) 

n—1 

= 2(xy)^e-f (^)+n,))^ P (^ ) ( ^ (vS) 

3=0 
n-1 

= 2(xy)^e^ v ^ +v ^ J2Pj N) ( x )pf ] (y) 

3=0 

= 2(xy)*K n;N (x,y). 
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This proves (12.31) . 

If a > and if we use (12.91) we obtain (12.41) in a similar way. □ 

We will also connect the density ipy of the equilibrium measure for V 
with the density ipw of the equilibrium measure for W. The equilibrium 
density ipy minimizes the energy (11.51) . It is characterized as the unique 
probability density function on [0, oo) with the property that for some 
constant Ey, 



log \x — y\ ipviy) dy — V(x) + ly = for x G supp r/v 5 

(2.11) 

2 / log \x - y\ ifo-(y) dy - V{x) + ly < for x G [0, oo). (2.12) 



The density ipw is the limiting mean eigenvalue density for the matrix 
ensemble ( 11.21) as n, N — > oo in such a way that n/N — > 1. Then %pw 
minimizes the weighted energy 



I w ^) = J J hg— — — ip(x)ip(y)dxdy + J W(x)ip(x) dx (2.13) 

among all probability density functions ip on M.. The density ipw is 
uniquely characterized by the following variational conditions for some 
constant £w- 

2 J log \x — y\ 4>w(y) dy — W(x) + i w = for x G supp ip w , 

(2.14) 



2 J log \x - y\ ifar(y) dy - W(x) + £ w < for iGl (2.15) 

Lemma 2.2. We have that 

ipw( x ) — klV'v'l^ 2 ) f or ^ G M. 

Proof. We first show that \x\ipy(x 2 ) is indeed a probability density. 
This follows from the positivity of ipy and the fact that 

roc roc 
\x\vjjy(x 2 )dx = 2 / Xll)y{x 2 )dx = I llly{s)ds = 1. 

> Jo Jo 

To prove part (i) of the lemma, it is now sufficient to show that 
\x\ipy(x 2 ) satisfies the variational conditions (I2.14p and (I2.15p . We 
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have 

\og\x-y\ \y\ip v {y 2 ) dy 



\og\x - y\yi) V {y 2 ) dy + \ log \x + y\ ytpv(y 2 ) dy 
\og\x 2 - y 2 \yi> v (y 2 )dy 



1 f°° 

- / log \x 2 — s\ ipv(s) ds. 

2 Jo 



This relation implies (l2~14j) and fl2~T5l by virtue of (l2~Tll and (l2~T2l) . 
Since these variational conditions uniquely characterize the lemma 
follows. □ 

In exactly the same way we find the following relation between the 
equilibrium densities of the sets Sy and Sw- Thus ujy is the unique 
probability density on Sy that minimizes the (unweighted) energy 

log io(x)uj{y) dxdy. 

\x - y\ 

Lemma 2.3. Let uiy be the equilibrium density of the set Sy and let 
u)w be the equilibrium density of Sy/ = {x \ x 2 G Sy}. Then 

co w (x) = \x\uj v (x 2 ) for x E S w . (2-16) 



3. Proof of Theorem 11.21 

The proof of Theorem 11.21 is based on the following result for the 
matrix ensemble (11.21) for the case that ipw vanishes quadratically at 0. 
There is a one-parameter family of limiting kernels that are connected 
with the Hastings-McLeod solution of the Painleve II equation q" = 
sq + 2g 3 — a. We denote these kernels by W^ lt,u (u, v ; s) and we refer 
to (I4.4p below for the precise form of K" lt,n (M, v ; s). The superscript 
II refers to the classification of singular points for the matrix ensemble 
(11.21) see [121 E] according to which there are three cases 

Singular case I: : Equality holds in (12.151) for some x G M. \ 
supp ipw Any such point x is then a singular exterior point in 
the external field W. 

Singular case II: : ipw vanishes at an interior point of its sup- 
port. Such a point is a singular interior point in the external 
field W. 
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Singular case III: : ipw vanishes at an edge point to higher 
order than a square root. Such a point is a singular edge point 
in the external field W. 

The following result was established in [Tj. 

Theorem 3.1. ([3 Theorem 1.2]) Let W be real analytic on R such 
that (II. 3p holds. Suppose thatipw{fy = V'w(O) = and ipwify > 0, and 
assume that there no other singular points besides 0. Define constants 

7T 7T 

ci,w = jrV'w(O) > anrf c 2,^ = -jj^-ujw(0) > (3.1) 

C 1,W 

where ujw is the density of the equilibrium measure of the support of 
ipw. Let n, N — ► oo such that lim n 2 ^ 3 (n/N — 1) = L E R exists and 

n,N—*oo 

put c = c^w/ anc? s = c 2t yyL. Then 

l r T^ K ffi(T^>T^) = ^( x >v> 8 ) (3 - 2) 

n,N-*oo {en) 1 ' 6 ' yicn) 1 ' 6 (cn) 1 > a J 
uniformly for x, y in compact subsets o/R \ {0}. 

Now we are ready for the proof of Theorem 11.21 with constants 

ci v = 77 lim \x~ 1/2 ip v (x)] and c 2 y = 4tt hm \«Jx uj v (x)] , 

C 1,V 

(3.3) 

where uiy is the density of the equilibrium measure of the support Sy 

Of Ipy. 

Proof of Theorem I j.£ We define 

W(a;) = W(x 2 ), xeR. 

Then by Lemma 12.21 we have that ipw{x) = \x\tpy(x 2 ), which implies 
by the assumptions on ipy near 0, that 

<M0) = ^(0) = 0, ^(0) = 2 lim [x-^tPvix)] > 0. (3.4) 

It follows from Lemma 12.21 that there are no other singular points be- 
sides in the external field W (in fact, the absence of type I singular 
points follows from the proof of the lemma). The conditions of Theo- 
rem [3J] are therefore satisfied. From ( 13. 31) . ( 13. II) and Lemmas 12.21 and 
it follows that 



cxy = 2c 1>w , and c 2 y = 2 2/3 



C2,W- 
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Then if n, N — > oo such that f ll.lOp and if c = C\y and s = CiyL y we 
get from from Theorem 13.11 that 

li m 1 K (2a+l,W) ( X y \ _ K crit,II f x 

and if a > we also have 

£ m ? u73 g 2n?2Jv' W) ( 7 ^1/3 W TTH ) = V, g )> ( 3 - 6 ) 

Then straightforward calculations based on (12.31) and (13.51) show that 
the limit (II. 9p exists with 

Kf /har V, y; s) = i(*y)- 1/4 (K^f ^ «) + ^f(v^, *)) > 

(3.7) 

while if a > 0, then from (12. 4p and (13.61) we also get 

C ft/har V, y; s) = \{xy)- llA ^ ») - d n (v^, - «)) • 

(3.8) 

□ 



4. Proof of Theorem 11.31 

To prove Theorem ll.3[ we recall the definition of the kernel IK" 1 *' 11 
from [7J. It uses the Hastings-McLeod solution q = q a of the Painleve 
II equation 

q" = sq + 2q 3 - a. (4.1) 
Given a > —1/2 and s G E, we abbreviate g = g Q (s) and r = g^(s), 
and we let \ ^"' l ^f'' \ \ be the unique solution of the equation 



d_ ( ® a ,i\ _ /-4^ 2 - i(s + 2q 2 ) Azq + 2ir + a/z\ ( $ Qil 

tfe V $ «:2 j " V AZ( 1 ~ ^ + Aiz2 + ^ + 2( f) J V $ 



Q,2 



(4.2) 



with asymptotics 



uniformly as z — > oo in the sector e < arg^ < it — e for any e > 0. 
Since q a has no real poles for a > —1/2, $ a i and $ Qj 2 are well-defined 
for real s. The functions <& a> i and $ Qj 2 extend to analytic functions on 
C \ (— ioo, 0], with branch points in 0. We denote these extensions also 
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by <& a ,i and <& a ,2- In the kernel K™*' 11 we need their values on the real 
line. Indeed, we have for real x, y, and s, 

lBr«*it,n/_ „. x _ Ania\ Sf rn(x)+s R n(y)] ®a,l{ X \ ^^{V^ S ) ~ s)^ a ,2{x] s) 

a V' y > s >- 2m{x-y) 

(4.4) 

We first rewrite IK™*' 11 in a 'real' form, using the functions (see also 
i), 

F a , 2 (x;s)J 2\-i i) \$ a , 2 (x;s)J ' X1J 

(4.5) 

Then we get from (I4.4p that 

^crit.ii/ „. x _ F a>1 (x; s)F ai2 (y; s) - F a>1 (y; s)F a , 2 (x; s) 
a ^ y > S >- 7r(x-y) 1 j 

and from (14.21) we get the differential equation 

A. ( Fa A - ( 4:xq + a/x Ax 2 + s + 2q 2 + 2r\ f F a \ 

dx~ \F^ 2 ) ~ \~4x 2 -s-2q 2 + 2r -Axq - a/x J \F a>2 

(4.7) 

The following properties are crucial. 

Lemma 4.1. Both F a>1 (x; s) and F a2 (x; s) are real for real x and s. 
In addition, F a ^ is an even function and F aj2 is an odd function, and 
they have the following asymptotic behaviors 

F a ,i(x; s) = cos (^x 3 + sx - + 0( x ~ l )i (4-8) 

F a>2 {x; s) = - sin ^x 3 + sx - T^aJ + ^ix" 1 ), (4.9) 
asi-> +oo. 

Proof. The proof is based on the symmetries of the Riemann-Hilbert 

problem that characterizes ( 

Let T = (J . Tj be the contour consisting of four straight rays oriented 
to infinity, where 

7T 57T 57T 

T 1 :&Ygz=-, r 2 :argz = — , T 3 : argz = — — , T 4 : arg^ 

DO 

see Figure |2j The contour T divides the complex plane into four sectors 
Sj, j = 1, 2, 3, 4 as shown in Figure [2j 

For a > —1/2 and s 6 1 we look for ty a (z;s) that satisfies the 
following Riemann-Hilbert problem. 
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r 3 



r 4 



Figure 2. The contour r consisting of four straight rays 
oriented to infinity. 



(a) ^ a : C \ T -> C 2x2 is analytic in C \ T. 

(b) ty a satisfies the following jump relations on T \ {0}, 



V a , + (z; s) = * Q _(z; s) ^-"L* , 

S) = * Q _(z; S) (^_^ ia ^ , 



1 



for z G I\, 
for z G T 2 , 

for z e r 3 , 

for z G T 4 . 



(4.10) 
(4.11) 
(4.12) 
(4.13) 



(c) \I> Q has the following behavior at infinity, 

V a (z] s) = (/ + 0(l/z))e- l( l z3+sz ^, asz^oo. 

Here c 3 = ( o -i ) denotes the third Pauli matrix. 

(d) ty a has the following behavior near the origin. If a < 0, 



(4.14) 





A 


a 


z\ 


c 


z 


a 


z\ 



(4.15) 
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and if a > 0, 



'o 


(\z\ 






\\A 




f\z\ 


< o 




\\A 


° 


f\z\ 







\z\ 
\z\ 



~ a \z\ a 



as z — > 0, z G Si U S3, 
as z — > 0, 2 G 5*2 , 
as 2 — > 0, 2; G £4. 



(4.16) 



It was shown in [7J that the Riemann-Hilbert problem for ty a has a 
unique solution and that 



^ a (z;s) 



VJz; s) 



^ a (z; s) 



^a(z; s) 



VJz; s) 



1 

-ma 







1 





for z E Si 

for z G S2, 
for z G S 3 , 

for z G S4, Rez > 0, 

for z G S4, Rez < 0. 



(4.17) 



It is easy to check that ai^f a (z; s)a\ and ai^? a (— z; s)o~i with a% = 
{il) also satisfy the conditions of the Riemann-Hilbert problem for 
\I/ a . Therefore by the uniqueness of the solution we have 

^ a (z; s) = cr^aiz; s)a x = <7i* tt (-z; s)^. (4.18) 

From the first equality in (14.181) it follows (as noted in Remark 2 of [7]) 
that 



e^ iassn{x) $ a:2 (x; s) = e^ iassnix) ® aA (x; s), for x G M\{0} and s G I 

(4.19) 

which shows that F a>1 (x; s) and F a)2 (x; s) are real, since we obtain from 
(g3D and (14TT91 that 

F aA (x; S ) = Re[e^ ias ^$ ail (x; S ) 
F a , 2 (x; s) = Im \e^ iassn{x) $ a ,i(x; s) 



(4.20) 
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From (14.141) . (14. 17ft . and (14.201) we obtain the asymptotic behaviors 
(Oil and 

The symmetry ^ a (z] s) = ai^f a (—z; s)ai yields by (I4.17P that 

e^ iasgn{x) ^ at2 (x; s) = e^ iassa (- x) <$> a>1 {-x; s), for x E R \ {0}, 

which in view of ( 14. 19j) and (14.201) shows that F a ^ is even and F a 2 is 
odd. ' □ 

Now we can give the proof of Theorem 11.31 

Proof of Theorem \1.3\ . From (13.71) and (14. 6p . we get 



jj^-soft/hard , 



x,y\s) 



1 , 1/4 f F a+hl (^)F a+h2 (^y) - F a+ll {^y)F a+h2 {^) 

2™ ^ AV^-Vy) 



7r(v^ + Vv) 



(4.21) 



Putting the fractions onto a common denominator, and using the even/odd 
properties of Lemma 14.11 we find 

Kr /hard (x )2 /; S ) 

1 ( W/4 ^a + il(v^)^a + i2(y^) -^ + l,l(v^)^ a+ l, 2 (v^ )' 

— \xy) I 

7r \ x — y 

(4.22) 

which leads to (11.111) with 

f a (x; s) = x-WF^i^VZ; s), (4.23) 

and 

g a (x;s) =x 1/4 F a+h2 (^;s). (4.24) 

The system of differential equations (11. 12ft follows from (14.71) . (14. 23ft . 
and (I4.24p and the asymptotic behavior (I1.14I) - (I1.15I) follows from (14.81) - 
(HD and 

This completes the proof of Theorem 11.31 □ 
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